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ABSTRACT 


This report is the User's Manual for the current CONIC 
computer program, written for GDB under Task. 427. This document 
revises and supersedes TRW Note 70-FMT-835 , CONIC Program User's 
Manual, May 27, 1970. The program has been written in FORTRAN 
for use on the TRW Time Share system. The CONIC program is 
used to calculate the relative motion of two bodies in any conic- 
section orbit (circular, elliptical, parabolic, or hyperbolic). 

This manual describes the basic equations used for trajectory 
propagation and defines the input and output parameters of the 
program. Also included are sample runs with corresponding figures 
and a detailed flow chart of the program. 
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1.0 INTRODUCTION 


The CONIC Program is designed for the calculation of the relative 
trajectories of two vehicles in orbit, where the orbits may be circular, 
elliptical, parabolic, or hyperbolic. The inertial trajectories of the 
two vehicles are calculated independently as conic sections, with no 
limitation on the type of orbit. 

The only inputs: required by the program are Cl) specification of the 
gravitational constant of the attracting body, (2) the initial positions 
and velocities of the primary and secondary vehicles, and (3) a list of 
the maneuvers to be performed by the secondary body. The maneuvers are 
specified by phase, where each phase consists of an impulsive maneuver of 
specific delta velocity and direction, followed by a nonpropulsive coast 
of specified duration. A maximum of ten such maneuver phases may be 
defined for a single run. 

The output of the program is the relative position and velocity of 
the secondary vehicle relative to the primary vehicle, printed at equal 
intervals during the coasting flight. The number of intervals per phase 
is controlled by program input. 

The program is ideally suited to the calculation of relative motion 
in a highly elliptical orbit such as a semi -synchronous earth orbit, or 
in a hyperbolic orbit such as a lunar flyby or an approach to earth entry 
It is also very accurate in computing the relative motion for complex man 
euvers in circular orbit, such as a rendezvous-rescue sequence. It is 
limited to impulsive maneuvers, and should not be used for long-term low- 
level thrusts such as continuous venting or atmospheric drag. 
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2.0 BASIC EQUATIONS 


The basic equation by which the trajectories of the bodies are 
calculated is the Kepler equation for nonpropul sive motion in a conic section 
(ellipse, parabola, or hyperbola), expressed in terms of universal variables. 
Use of this equation allows direct solution of the problem of predicting the 
position and velocity of a body after a given time interval n t Tt , given its 
initial position and velocity. Kepler T s equation may be written as: 



— X 2 C(a X 2 ) + (1 
o 


r a )x 3 S(cc x 2 ) + r x 
O O 0 o 


where the functions S(£) and C(t>) are defined as: 
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and 

|i is the gravitational constant 

t is the time interval of propagation 

r is the initial position vector of the body 

o r 

v is the initial velocity vector of the body 
o 

is the magnitude of r^ 

v is the magnitude of v 

o o 

x is the normalized anomaly 

An examination of Kepler r s equation shows that it is an explicit expression 
for the known time of propagation 11 1 11 in terms of the unknown variable "x", 
To evaluate x for a given value of t requires an inversion of the function 
t = t(x) given by Equation 1. This is done in a straightforward manner by 
employing Newton r s method, as follows: given an estimate = x n ^ t n ) t * ie 

required value of x “ x(t), an improved estimate x ^ is obtained from: 



where the derivative is given by; 



a Q X 3 S(a o x 2 )] 
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( 6 ) 


The process is repeated until a sufficiently accurate value of x is obtained. 
Once x has been determined, the position and velocity vectors of the body at 
time t are given by: 
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where r is the magnitude of r(t). 


( 8 ) 
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The foregoing method for propagating the state vector of a body is 
applied independently to the primary and secondary bodies to obtain their 
inertial state vectors at the required times. The relative motion parameters 
are then obtained by subtracting the state vector of the primary body from 
that of the secondary body, and expressing the resultant relative vector in 
a local •vertical -local horizontal system based on the primary vehicle* The 
process is repeated for each printout in a given coast phase. At the 
beginning of a new phase, the delta velocity vector of the impulsive maneuver 
is converted to inertial coordinates and added to the inertial velocity 
of the secondary vehicle. The resulting state vectors of the primary and 
secondary bodies are then propagated for the duration of the phase. 
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3.0 INPUT PARAMETERS 


A total of seven input options have been included in CONIC to allow 
for maximum flexibility in providing input to the program. These options 
are as follows: 

Input option 1 (INOPTl) allows a choice of coordinate systems. If 
IN0PT1 is set equal to 1, the relative position and velocity of the 
secondary body is expressed in a local-vertical-local -horizontal Cartesian 
coordinate system centered on the primary body. If INOPTl is set equal 
to 2, the relative position and velocity of the primary and secondary 
bodies are converted to state vectors in an inertial coordinate system. 

The second input option (INOPT) allows a choice of methods for specify- 
ing the pitch direction of maneuvers executed by the secondary body. If 
INOPT is set equal to 1, the pitch maneuvers are measured upward from the 
direction of the posigrade local horizontal. This method of specifying 
the pitch of a maneuver is used both for elliptical orbits and for hyperbolic 
orbits. If INOPT is set equal to 2, the pitch maneuvers are measured upward 
from the direction of the inertial velocity vector of the secondary body. 

This method of specifying the pitch of a maneuver is most useful for hyperbolic 
orbits such as lunar flyby or earth entry orbits. 

Input option 2 (IN0PT2) allows a choice of output units. If IN0PT2 
is set equal to 1, the relative position of the secondary body will be 
typed out in nautical miles. If IN0PT2 is set equal to 2, it will be typed 
out in feet. 

Input option 3 (IN0PT3) is used to define the gravitational constant 
(EMU) and radius (R0) of the central gravitational body. If IN0PT3 is 
set at 1, an earth orbit is assumed, and if IN0PT3 is set at 2, a lunar 
orbit is assumed. Direct input of EMU and RO Is not required, because 
the necessary values are included in the program. 
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Input option L, (IITOFT4) is used to choose the method for initializing 
the position and velocity of the primary (or reference) body. Setting 
IN0PT4 equal to 1 allows for direct input of the inertial position (XOP) 
and velocity (VOP) vectors of the primary body in feet and feet per second. 
If IN0PT4 is set equal to 2, the primary body is assumed to be in a circular 
orbit, and the only required input parameter is the altitude of the circular 
orbit (H) in nautical miles. The state vector for the primary body is then 
calculated automatically according to the formulas: 

X0P(1) = 0.0 V0P(1) = SQRT(EMU/X0P(2)) 

XOP (2) = RO + H * 6076.1155 V0P(2) =0.0 

XOP (3) = 0.0 VOP (3) = 0.0 

Regardless of how the initial state vector for the primary body is input, 
it is used to calculate the type of orbit in which the primary body will 
move. If the orbit is hyperbolic or parabolic, a message is typed to 
identify the orbit type. If the orbit is elliptical, the period is 
computed and typed. The orbital period is needed when the coast times 
between maneuvers are known only in orbit fractions. 

Input option 5 (IN0PT5) allows a choice of methods for initializing 
the position (XOS) and velocity (VOS) of the secondary body. If IN0PT5 
is set equal to 1, the initial position and velocity of the secondary body 
are set equal to that of the primary body. If IN0PT5 is set equal to 2, 
the initial position and velocity of the secondary body may be specified 
relative to the primary body in a local-vertical-local-horizontal coordi- 
nate system based on the primary body. This coordinate system is defined 
with the "+X" axis in the local horizontal posigrade direction, the "+Y" 
axis in the local vertical upward direction, and the "Z" axis completing 
a right-handed system. If IN0PT5 is set equal to 3, the inertial state 
vector of the secondary body may be input directly. This option is 
available only if the primary body was also initialized by direct speci- 
fication of the inertial state vector (IN0PT4 = 1). 
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Input option 6 (IN0PT6) is encountered at the end of each case, and 
provides a choice of three ways to begin the next case. When the program 
is executed from statement 1, the first required input is IN0PT3, which 
allows the gravitational body to be specified as either the earth or the 
moon. The next required inputs are the initial state vectors of the two 
bodies. The last required inputs are the maneuvers to be executed by the 
secondary body. If IN0PT6 is set equal to 1, execution of the next case 
begins with the input of new initial state vectors and continues with the 
input of a new table of maneuvers. If IN0PT6 is set equal to 2, execution 
of the next case begins with the input of a new table of maneuvers, leaving 
the initial state vectors unchanged. If IN0PT6 is set equal to 3, the 
program branches to STOP, This option would be used when no more cases 
are to be executed, or when the next case must begin execution at statement 
1 in order to change the value of IN0PT3, which specifies the gravitational 
body. 


In addition to the inputs provided for by the various input options, 
the maneuvers to be executed by the primary or secondary body must also be 
supplied. Six parameters must be input for each phase, where each phase 
consists of an impulsive maneuver followed by coasting flight. The impulsive 
maneuver is executed either by the primary or secondary body (B) . The 
impulsive maneuver is defined by the pitch (P) and yaw (Y) directions of the 
maneuver and its magnitude (DV) . The pitch (P) is input in degrees, measured 
upward from either the local horizontal posigrade or the inertial velocity 
vector, depending on whether INOPT was set equal to 1 or 2, respectively. 

The yaw (Y) is also Input in degrees, measured positively to the right of 
the orbital plane (facing in the direction of the inertial velocity) . The 
magnitude of the impulsive maneuver (DV) is input in feet per second. The 
coasting flight following the impulsive maneuver is specified by the total 
coast time (CT) , input in seconds, and the number of points (NP) to be 
plotted per phase. The total input for each phase consists of the primary 
or secondary body (B) , pitch (P) , yaw (Y) , and delta velocity (DV) of the 
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impulsive maneuver, and the coast time (CT) and number of points to be 
plotted (NP) for the subsequent coasting flight. The inputs for pitch 
(P) , yaw (Y) , delta velocity (DV) , coast time (CT) , and number of points 
(NP) must be in decimal form. Any number of phases may be specified 
(NPHASE), up to a maximum of ten phases. The input for the specific 
number of maneuver phases must be in 13 (INTEGER) format, i.e., one 
maneuver phase is inputted as 001. 
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4 . 0 OUTPUT PARAMET ERS 


The output parameters for CONIC consist of the relative position and 
velocity of the secondary body 'expressed in a local-vertical-local-horizontal 
Cartesian coordinate system centered on the primary body. The definitions 
of the output parameters are as follows: 

T Time duration measured from the beginning of the first 

phase (seconds) 

RPDR Relative position downrange of the secondary body, positive 

in the direction of the inertial velocity of the primary 
body (nautical miles or feet) 

RPRA...... Relative position radial of the secondary body, positive 

upward (nautical miles or feet) 

RPCR Relative position crossrange of the secondary body, positive 

to the right of the orbital plane when facing in the direction 
of the inertial velocity (nautical miles or feet) 

RVDR Relative velocity downrange of the secondary body (feet per 

second) 

RVRA Relative velocity radially up of the secondary body (feet per 

second) 

RVCR.,.,.. Relative velocity crossrange of the secondary body (feet per 
second) 
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5.0 SAMPLE COMPUTER RUNS 

The following listings and figures show the input and output of several 
computer runs. The output data is plotted on Figures 1 and 2 for the 
corresponding inputs. The input parameters required by the program and 
supplied by the user are underlined; all of the remaining print-out is 
typed automatically by the program. The input options for Figure 1 are 
the same as those for Figure 2. The first four input options are set equal 
to 1. This sets the options so that a local-vertical-local-horizontal 
coordinate system is assumed, pitch maneuvers are measured upward from the 
local horizontal posigrade, the relative position is output in nautical 
miles, and so that an earth orbit is assumed. The primary body is 
initialized by specification of its circular orbit altitude. The secondary 
body is initialized with its relative position and velocity equal to that 
of the primary body. 

Immediately after initialization of the primary body, its type of orbit 
is calculated, and the orbit type is identified if it is hyperbolic or 
parabolic. The orbital period is also calculated and printed out. The 
period is needed because coast times between maneuvers must be input in 
seconds. This allows the specification of coast times equal to a period or 
even to a fraction of a period. All of these sample runs are executed for 
one maneuver phase. 

The maneuver table allows a detailed specification of the maneuvers to 
be performed either by the primary or secondary body. The maneuvers are 
specified by phases, where each phase consists of an impulsive maneuver, 
followed by a coasting flight during which the relative position and velocity 
of the specific body are printed out. The maneuvers in both Figures 1 and 2 
were performed on a coasting flight of 10,000.0 seconds with printout at 
10 equal intervals during the coast. 

The output blocks show the result of the input parameters. The first 
line of output shows the relative position and velocity of the secondary 
vehicle before the first maneuver. Since the secondary vehicle was ini- 
tialized with its position and velocity equal to that of the primary vehicle. 
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its relative position and velocity are zero* The second line of output 
(the first line under Phase Number 1) shows the relative position and 
velocity of the secondary vehicle immediately after the phase 1 maneuver. 

The remaining lines of output under phase 1 give the relative position and 
velocity of the secondary body for ten equal increments during a one-orbit 
coast. The listings of maneuvers I, Ha, and Ilia in Figure 1 were ran- 
domly selected as sample print-outs. 

Three sample posigrade separation maneuvers at two different delta 
velocities are plotted in Figure 1. All three maneuvers were executed in- 
plane (0 degree yaw) and at 3.3 feet per second and 7.0 feet per second, 
respectively. In maneuvers I and la, a 0 degree pitch is specified. A 
pitch of 45 degrees upward is specified in maneuvers II and Ila, A pitch 
of 90 degrees (radially upward) is specified in maneuvers III and Ilia. The 
relative motions of the last two maneuvers return to the original position 
at the end of one orbit resulting in football maneuvers. Figure 1 shows 
only the side views of the relative motion of the several maneuvers. 

Two sample posigrade separation maneuvers at different delta velocities 
are plotted in Figure 2. Both maneuvers were executed with 30 degrees pitch 
upward and 30 degrees yaw out of plane. Maneuver I was performed at 3.3 
feet per second and Maneuver II at 7.0 feet per second. Figure 2 shows both 
the side and top views of the relative motion of the two maneuvers. 
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Note that the relative position is still zero, but the relative velocity 
is now upward and out of plane, as specified by the input. The remaining 
lines of output under phase 1 give the relative position and velocity of 
the secondary body for ten equal increments during a one-orbit coast. The 
in-plane motion is a football maneuver, returning to the original position 
at the end of one orbit. 

In the second run, all of the input options were set equal 2, so that 
a lunar orbit is assumed, the primary body is initialized by specification 
of its circular orbit altitude, the secondary body is initialized by specifying 
its position and velocity relative to the primary body, the pitch of maneuvers 
is measured upward from the inertial velocity vector of the secondary body, 
and so that relative position is output in feet. The maneuver specified is 
10 feet per second radially up (football maneuver), with output at ten points 
equally spaced over one orbit. 
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INPUT OPTIONS for FIGURES 1 and 2 
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FIGURE 1 - MANEUVER Ila 
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Ei Figure 2. Motion of a Secondary Body Relative to a Primary Body - 
Posigrade Separation Maneuvers (Side and Top Views) 
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6.0 DETAILED PROGRAM FLOW CHART 




CONT. ON PO 2 










RERD ( 5.2) I NOP T3 


IF( I N0PT3- 1 ) 8 . 8 . 9 























com. ON PG 4 


6-3 














ARE = SQRT ( XOPC L )**2+X0P(2)**2+X0P(3)**2) 


VEE = SQRT(V0P( l )w*2+VOP(2)*»2+V0P(3)**2) 
ENE = 0.5*VEE**2 - EMU/ARE 


C0NT. ON PG 5 


6-4 




















I F( IF [ X ( ENE) ) 26.24,22 


22 V 





PERIOD = 2 -0*3- l415926535*SQRT(fl«w3/EMU 
WR I TE( 6 . 27 ) PERIOD 




















CONT. ON PG 7 


6-6 











FORMAT! * IN0PT5C INITIALIZATION OPTION FOR SECONOARY BOOT)*) 


IF ( IN0PT4-1 J33.33.35 


32 

WRITE 


FORMAT ( *1 *.SAME AS PRIMARY .2 ! INPUT RELATIVE VECTOR .* 
» *3 ; INPUT STATE VECTOR*) 


READ! 5.13) IN0PT5 


IF ( IN0PT5-2J37.42.38 


35 

I WRITE (6. 367 


FORMAT! *1 ‘.SAME AS PR IMARY .2 : INPUT RELATIVE VECTOR*) 


REA0!5. 13 ) IN0PT5 


IF ( IN0PT5- 1 J37.37.42 


37 <7 

XOS! 1 ) = XOP( 1 ) 


CONT. ON PG 8 


6-7 



































CONT. ON PO 10 


6-9 
















54 



CONT. ON PO 11 


6-1 
















CONT. ON PG 12 


6-11 













CONT. ON PG 13 

6-12 












CONT. ON PG 14 


6-13 








V0SR ( 2 )=VOSR( l }*SIN(THETR)+VOSR(2)*C0S(TNETR) 
VOSR(n = RHO 


n 


B2 


-cH 


91 


XOSR( l )=XOSR(l )*OUT 


XOSR(2)=X0SR(2)*0UT 
XOSRt 3) = X0SR(3)*OUr 

WRITE(6.102)T.X0SR( 1) .XOSRt 2 1 .XOSR t 3 ) . VOSRt i) . V0SR( 2 ) .VOSRt 3 ) 


102 

F0RHRnF7.0.3Fl0-l . F 1 1 .3.2F1Q.3) 

z s t 

< DO 121 1 = 1 ,NPHRSE > - - 


TPHRSE=0.0 
EXPE=0 .0 
EXES=0 .0 
MRlTEtS. 103)1 


l> 19 


103 <7 



C0NT . ON PG 15 


6-14 










C0NT. ON PG 16 


6-15 











CONT. ON PG 17 


6-16 











^ B6 

115 

WR[TE(6,102)T,X0SR( 1 ) ,X0SR ( 2 ) .XOSR ( 3 ) . VOSR ( 1 ) ,VOSR( 2 ) .VOSR( 3 ) 


M 


DT=CT( I ) /NP ( I ) 
K=NP ( I ) 


00 120 J=1 .K > - - - [> 18 


T=T+OT 

TPNflSE=TPHflSE+DT 




CONT. ON PG 18 


6-17 















CONT. ON PG 19 


6-18 













14 


> i 


VS( 3 ) = ' 

VES( 3 ) 

b 

WRITECi 

B.200) 


200 v 

FORMAT (*IN0PT6( 1 :NEW VECTORS AND NEW MANEUVERS , 2 :NEW* 
,* MANEUVERS ONLY >3:ST0R )») 


5 


READ ( 5.2C 

11 ) IN0PT6 

on 1 r 

T 


201 <7 



END 


6-19 










CONT. ON PCr 2 


6-20 











END 


FUNCTION SINHCX) 

E=2 .718281828459045 
SINH=0.5*(E**X-E**l-X ) ) 


RETURN 
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6-23 





















